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In this paper, we provide a simplified stellar structure model for ideal gases, in 
which the particles are only driven by gravity. According to the model, the structural 
information of the star can be roughly solved by the total mass and radius of a star. 
To get more accurate results, the model should be modified by introducing other 
interaction among particles and rotation of the star. 
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I. SOME PHENOMENA INSIDE A STAR 

The spherical symmetric metric for a static star is described by Schwarzschild metric 

g^u = diag (b(r), -a(r), -r 2 , -r 2 sin 2 9) . (1.1) 

For the energy momentum tensor of perfect fluid 

Tfo, = (p + P)UJJ V - JV = diag (bp, aP, r 2 P, r 2 sin 2 6P) , (1.2) 

where p(r),P(r) are proper mass energy density and pressure, = (^,0,0,0), we have 
the independent equations as follow 

(£)' = 1 - SirGpr 2 , (G 00 = -87rGT 00 ), (1.3) 

h - = — - + 8nGPar, (G u = -8tiGT u ), (1.4) 
o r 

P ' = -(P + P )& P* = 0). (1.5) 

To determine the stellar structure of an irrotational star, we solve these equations. How- 
ever (1.3)-(1.4) is not a closed system, the solution depends on an equation of state(EOS) 
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P = P(p). For polytropes, we take P = P P 7 [1]- For the compact stars, there are a lot of 
EOS derived from particle models [2] -[7], which provide the structural information and pa- 
rameters such as the maximum mass for neutron stars. However, in these EOS the boosting 
effects of the gravity on particles seem to be overlooked or adopted unconsciously. As shown 
in [8], the static equilibrium equation (1.5) is insufficient to describe such dynamical effects 
on fluid, and a Gibbs' type law should be included. In this paper, according to this Gibbs' 
type law, we derive the stellar structure model, which shows that the boosting effects of 
gravity play an important role to the stellar structure. 

Before expanding the model, we make a few simple calculations and examine the behavior 
of the metric and particles inside a star to get some intuition. The first phenomenon is that, 
the temporal singularity and spatial singularity occur at different time and place if the 
spacetime becomes singular, and the temporal one seems to occur firstly. 

Denoting the mass distribution by 

nr 

M(r) = AnG / pr 2 dr, K = 2M(r). (1.6) 
Jo 

Then by (1.3) we have solution 







(1- 


r J 
r ) 



;i.7) 



, if r < R, 
if r > R, 

where R is the radius of the star, and the Schwarzschild radius becomes 

r R 

R s = 2M(R) = K(R) = 8nG / pr 2 dr. (1.8) 

Jo 

For any normal star with R > R s . From the above solution we learn a(r) is a continuous 
function and 

a > 1, a(0) = lim a = 1, a max = a(r m ), (0 < r m < R). (1.9) 

r-+oo 

So the spatial singularity a — > oo does not appear at the center of the star when the 
singularity begins to form. 

On the other hand, by (1.4) and P > 0, we find b'(r) is a continuous function satisfying 

&'(()) = 0, b'(r) > 0, (Vr > 0), 6 = 1-^, (r>R). (1.10) 

(1.9) shows b(r) is a monotone increasing function of r with smoothness at least C 1 ([0, oo)). 
Consequently, the temporal singularity 6^0 should take place at the center. 
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FIG. 1: The trends of (a, b) as R — > R s , which show the spatial singularity does not occur at the 
center, but temporal singularity may occur at the center before a — > oo 

The trends of a(r), 6(r) are shown in FIG.l, where we take the Schwarzschild radius 
R s = 1 as length unit. From FIG.l and some simplified calculations, we find 6(0) — > seems 
to appear before the space becomes singular a — > oo. 

The second phenomenon is that, the particles near the center of the star are unbalanced, 
and violent explosion takes place inside the star before the temporal singularity occurs. 
When 6(0) -> 0, by (1.10) we have 

b^b r a , (a>l), if r < R. (1.11) 

Substituting it into (1.5), we find 

OL 

-P'^(p + P)-^+oo, (r-0). (1.12) 

According to fluid mechanics, —d r P corresponds to the radial boosting force, so (1.12) means 
violent explosion. 

More clearly, we examine the motion of a particle inside the star. Solving the geodesic 
in the orthogonal subspace (t,r,9), we get [9, 10] 

t= cTV 9 = ^> ^ = °' r =r{w*-^- 1 )' (113) 
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where Ci, C 2 are constants. The normal velocity of the particle is given by 




adr 2 



'K i+ ^)' 




r 2 d6 2 _ C 2 C 2 b 




(1.14) 



Mt 2 



Mt 2 



The sum of the speeds provides an equality similar to the energy conservation law 



(1.15) holds for all particles with v v ^ due to the symmetry of the spacetime. 

From (1.15) we learn v — > 1 when b — > 0, this means all particles escape at light velocity 
when the temporal singularity occurs. So instead of a final collapse, the fate of a star with 
heavy mass may be explosion and disintegration. The gravity of a star drives the inside 
particles to move rapidly and leads to high temperature. How the particles to react to the 
collapse of a star needs further research with dynamical models. A heuristic computation 
for axisymmetrical collapse is presented in [11], which reveals that the fate of a collapsing 
star sensitively depends on the parameters in the EOS. 



In this paper, we simply take the star as a ball of ideal gases, which satisfies the following 
assumptions: 

(Al) All particles are classical ones only driven by the gravity, namely, they are charac- 
terized by 4-vector momentum p 1 ^ and move along geodesic. 

(A2) The collisions among particles are elastic, and then they can be ignored in statistical 
sense [9, 10]. 

(A3) The nuclear reaction and radiation are stable and slowly varying process in a normal 
star, in contrast with the mass energy, the energy related to this process is small noise, so 
we treat all photons as particles and omit the process of its generation and radiation. 

These are some usual assumptions suitable for fluid stars. However, together with (1.3)- 
(1.4), they are enough to give us a simplified self consistent stellar structure theory. In [8], 
we derived the EOS for such system as follows 






II. THE EQUATIONS FOR STELLAR STRUCTURE 



p 



3 k 
Af [J (J + 2<j)]2 , , 

mc 2 

x(fhc 2 + hT) = Q [J(J + 2a)] 3 * (l + lA c 2 , 




(2.2) 



(2.1) 



P 



(2.3) 
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where Af is the number density of particles, Ao = A/o(m, a) is related to property of the 
particles but independent of J, c = 2.99 x 10 8 m/s the light velocity, <r=| a factor reflecting 
the energy distribution function, m the mean static mass of all particles, J dimensionless 
temperature, which is used as independent variable, (p, P) are usual energy density and 
pressure, g a mass density defined by 

g = A/" m. (2.4) 

By (2.2) and (2.3), we get the polytropic index 7 is not a constant for large range of T 
satisfying 1 < 7 < |, and the velocity of sound 

[dP _V3 { c 2 J{2a + J) (5a 2 + 8a J + 4 J 2 ) \ * ^3 
Gsound " C y dp ~ 3 V(a + J)2[2<7 + (2 + 5a)J + 4J2]yl < 3°' ( j 

which shows that the EOS is regular and the causality condition holds. 
Equation (1.5) can be rewritten as 

db 2b dP 



dJ~ p + PdJ' 
Substituting (2.2) and (2.3) into (2.6), we solve 

A(R - R s )(a + J) 2 



R[2a+ (2 + 5a)J + AJ 2 } 2 ' 



(2.6) 



(2.7) 



where R is the radius of the star, and R s = 2M tot the Schwarzschild radius. Substituting 
(1.7) and (2.2)-(2.7) into (1.3) and (1.4), instead of the Tolman-Oppenheimer-Volkoff equa- 
tion and Lane-Emden equation [1], we get the following dimensionless equations for stellar 
structure, 

n\r) = Q (2 + 3J)[J(J + 2(j)]i, (2.8) 

J/ W = 2^((0 2 [ J ( J + 2 ^ f + > + J ))' ^ 
where x is a constant length scale, 

X = c(AnGg)-^ = c(AnGM m)-^, (2.10) 

_ 2a + (2 + 5a)J + 4J 2 
^ = 5a 2 + 8aJ + A.P ' (C ~ ^ (2 ' 11} 

The exact solution to (2.8) and (2.9) seems not easy to be obtained. However, they are 

dimensionless equations convenient for numeric resolution. If we take x = 1 as the unit of 

length, the solution can be uniquely determined by the following boundary conditions 

K(0) = 0, J(0) = J > 0, J{r) = 0, (Vr > R). (2.12) 
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Dimensionless Mass vs. Temperature 

T 



Dimensionless Mass vs. Radius 




(Length unit: x) 
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FIG. 2: Relation between mass and cen- 
tral temperature similar to H-R dia- 
gram, where it is concentrated by the 
scale x an d J- 



FIG. 3: Relation between mass and Ra- 
dius. All structural information of a 
star are determined by a given radii pair 
(R s , R). 



Adjusting J , we get different solution. Some solutions are displayed in FIG.2-FIG.5. Usually 
we can easily measure the radius R and the total mass or equivalent Schwarzschild radius 
R s of a star. Then the other structural parameters can be determined by this radii pair 
(R s , R). In what follows, we show how to use FIG. 2 and FIG. 3 to solve practical problems. 
For the sun, we have the radii pair as 

R s = 2.96 X 10 3 m, R Q = 6.96 x I0 8 m. (2.13) 

By R s /R & = 4.25 x I0~ 6 , according to the relation shown in FIG. 3, we can solve the 
intersection A and get the dimensionless radii pair 

X = \og(R/ X ) = 2.320, Y = \og(R s / X ) = -3.052. (2.14) 
Consequently, we have 

x = 10~ 2 - 32 R = 3.335 x 10 6 m. (2.15) 

By Y we get intersection B in FIG. 2, and then get the central dimensionless temperature 
Jo = 1.145 x 10~ 6 for the sun. Taking it as initial value we can solve (2.8) and (2.9), and 
then get detailed structural information for the sun. 
By (2.10) and (2.15), we get 

c 2 

g = = 9.640 x 10 12 (kg/m 3 ). (2.16) 
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Influence of o on Mass M 



Dimension less functions of Neutron Star 
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bution on solutions. The results are not compactest stars. The trends are typi- 



sensitive to a. 



cal for all stars. 



Then by (2.2) and (2.3) we solve the mass density and pressure at the center 

p(0) = q[J (J + 2<rp (l + | Jo) = 8.45 x 10 3 (kg/m 3 ), 

P(0) = ^g[J {J + 2(j)]i((T + JorV = 8.70 x 10 8 (MPa). 

The temperature depends on the mean mass fh. By (2.1), we have 

T c = mc 2 J /k = n p {m p c 2 J /k) = 1.247 x 10 7 n p (K), 

where m p = 1.673 x 10 _27 kg is the static mass of proton, n p is the equivalent proton number 
for the particles. 

If the ionization in the sun is about H + + iV + + 2e~, then we have 



(2.17) 
(2.18) 

(2.19) 



n p = (70% x 1 + 30% x 14)/4 = 1.23, 

fh = n p m p = 2.06 x 10~ 27 (kg), 

T c = 1.247 x 10 7 n p = 1.53 x 10 7 (K), 

W = c 2 /(4vrGmx 2 ) = 4.68 x 10 39 (m^ 3 ). 



(2.20) 
(2.21) 
(2.22) 
(2.23) 



The compactest stars(with the maximum R s /R) correspond to the points M, N in FIG. 3 
and FIG. 2. The radii pair is R s : R = 1.59 : 5.77 and the central temperature Jq = 0.30. 
For a compact star with the solar mass M Q , we get the length unit and radius as 



X = R a / 1.59 = 1.86 km, R = 5.77x = 10.7 km. 



(2.24) 
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Along the above procedure, we solve 

p(0) = 8.50 x 10 18 kg/m 3 , (2.25) 
P(0) = 1.24 x 10 29 MPa, (2.26) 
T(0) = 3.27 x 10 12 n p K. (2.27) 

These are typical data for a neutron star[l]-[7]. The metric functions (a, b) and the mass, 
temperature distributions (7Z, J) for this star are displayed in FIG. 5. 

III. DISCUSSION AND CONCLUSION 

The above numerical results show that, although the orbits of the particles are not simple 
ellipses due to collisions, the influence of the gravity still exists and leads to high temperature 
inside the star. The main part of the EOS may be related to gravity. 

In contrast with (2.17) and (2.18), we find the central density and pressure in the sun are 
about one order of magnitude less than the current data 

p(0) = 1.6 x 10 5 kg/m 3 , P(0) = 2.5 x 10 10 MPa. (3.1) 

This difference is an unsettled serious problem, which seems to be caused by the dynamical 
effect of gravity. 

FIG. 4 shows that the solutions are not very sensitive to the concrete energy distribu- 
tion or a, so one needs not to solve specific problems via calculating complex distribution 
functions [8]. 

For photons, by the Stefan-Boltzmann's law 



we get 




How to determine the concrete function Af (m, a) is an interesting problem. 

The dimensionless equations (2.8) and (2.9) simplify the relations between parameters. 
This function is similar to that of the similarity theory [12]. However in these equations the 
information of the interaction among particles is ignored, so it can not provide the critical 
data such as the maximum density, the largest mass. To get such data the potentials and 
interactive fields should be introduced to the energy momentum tensor [8]. 
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